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I 


Let Aji-1, 2,3, 4) and B,; be eight integers such that 


SAY =SB" (n=1,3,5). 


ist i=1 
Then, if we put 
1=x—y, As=yt+w, As=z—v, Ag=v+u4, 
B,=x+y, B=w—y, B=z+v, Bg, =u—v, (2) 
we have to determine x, y, z, u, v and w soas to satisfy the Biophantine system 
v (u® — x?) = yp (x? — w?) (2v)? + (u+ z)? + (u — 2)? 
= (2y)? + (x + w)? + (x — w)?. (3) 
Now, by means of Goldbach’s identity 
p* + g* — Gr — p — 9g)? = Cr — p)* + (2r — @)* 
+ (r — p — 4)’, 
we see that Eq. (3,) is satisfied if we have 
v=tp,ut+z=3r—p—q,u—z=q 
y=4(2r—p), x+w=2r—q, x-w=r—p—g. 
Hence we get 
u = 3 (3r — p) 
z= x =4(3r— p — 29) 
w=4(r+p) 
and the first of Eq. (3) takes the form 
2r? — 3(p +9) r + (yp? + pg + gq’) = 0. 
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Therefore, 


p= 3D +9) + VP? + 10pq + @? 
= 4 


and, to exclude the case r = irrational, it is necessary and sufficient to take 
p? + 10pq + ge = 0. 
From this we find 
p = 2ab + 1067, q=a*® — b*, t= a* + 10ab + b* 
and then we have 
r, = a? + 4ab + 7b?, rg = 4 (a? — 2ab + 135°). 


Substituting the foregoing values in the Eqs. (5,), (54) and (6) we get 


xX, = Zz, = 4 (a + 10ab + 135?) Xy = Zz = 4(— a? — 10ab + 235°) 
y, = a? + 3ab + 25? Vo = 4 (a? — 4ab + 3b?) 
w, = $(a? + 6ab + 176?) We = 4 (a? + 2ab + 33b?) 
v, = ab+ 5b? Vo = ab + Sb? 
u, = 4 (3a? + 10ab + 115?) Us = 4 (3a? — 10ab + 195’). 





Thus, from (2) and (1), we see that either 
— qa? + 4ab + 9b", 3a? + 12ab + 21b*, a? + 8ab + 362, 
3a? + 12ab + 215? 
= 3a? + 1l6ab + 17b?, — a? + 13b2, a® + 12ab + 2352, 
3a* + 8ab + b*, (where n = 1, 3, 5), (7*) 


or 


— 3a* — 2ab + 176%, 3a? — 6ab + 39b?, — a? — 14ab + 35%, 
3a? — bab + 39b? 
= a? — 18ab + 29b?, — a® + 10ab + 27b?, 
— a® — 6ab + 43b*, 3a? — 14ab—b? (n=1, 3, 5). 
(8) 


* The separation of two sets of numbers by the symbol = denotes that they have the same 
sums of kth power for k = 1, 2,..n, or for noted values of x. 
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Hence we have a process for finding solutions of the system (1) and 
in the case: A,=A,, A, + A, = B, + By. 
Also, by means of (3), we find that either 
(2ab + 10b*) (2a? + 10ab +- 12b*) (a? — b?) | 
= (2a? + 6ab + 4b?) (a? + 8ab + 155?) | 





x (2ab — 26%) 
(2ab + 10b*)? + (2a? + 10ab + 12b?)? + (a? — b*)? - 
= (2a* + 6ab + 46%)? + (a? + 8ab + 15b*)? | 
+ (2ab — 25%)? 
or 
(2ab +- 10b?) (et a) (a? — b?) | 
= (a? + 4ab + 3b*) (— 2ab + 146?) 
- — a — Gab — =) | 
| 
(10) 


(Qab + 10b%)? + (< et ay + (a? — 69) 


= (a? — 4ab + 3b)? + (— 2ab + 1462)? | 


; é a — 6ab — 8 


~ 


It is easy to see that the equalities (9) express the condition that two 
rectangular parallelopipeds shall have integral edges, equal volumes and 
equal diagonals, if all the factors are positive. 


Though the linear substitutions a= a,— 5b, b=2b, transform, 
finally, the relation (7) into (8), however we consider them both, in order 
to find some other results. Also we consider (9) and (10). 

II 


From (4) we get the double relation 


Ps 3r —p—q= r—p, 2r—q, p+q-—r. 
Hence (9,) and (10,) become 
2ab + 1067, 2a? + 10ab + 12b7, a? — b? 
= 2a? + 6ab + 4b*, a® + 8ab + 1562, 2b? — 2ab 
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2 2 
a? — 4ab + 3b, — 2ab + 14b2, * = ~~ == 


2 2 
2 2b + 1059, ©— SOT ow, 


from which we find 


2 
2a?+10ab+12b?, a?—4ab+3b%, —2ab-+14b?, — —= 


= 2a? + 6ab + 4b%, a® + 8ab + 15b*, — 2ab + 262, 


a* — 10ab + 215? 
= ; 


On the other hand, from (9,) and (10,), we take 
(2a*-+-10ab+-12b%) (a®—4ab-+3b2) (—2ab-+14b4 (+606 + 5b) 





= (2a? + 6ab + 4b?) (a® + 8ab + 15b*) (— 2ab + 2b?) 


x (Se 4 0) 


Therefore we have a process for finding solutions of the system 


A,, As, As, A, = By, Ba, Bs, Bg 
A,A,AsA, = B,B.B,B,. 

Thus, 
4, 8, 30, 35 = 3, 14, 20, 40 
4.8.30.35 = 3.14.20.40 


and the polynomials 
x* — 77x® + 1862x? — 14680x + 33600, 
x4 + 77x + 1862x? — 16120x + 33600 
have respectively the roots (4, 8, 30, 35) and (3, 14, 20, 40). 


iit 


A two-parameter solution of the system 
XytXgtxXgt XH At wt wt % 


Nyt Xq=X%etXye=H 1 t+m~=— wt 7%, (11) 


XyXoXgX4 = 23292324 


may be found from (7). 
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In fact, for A, + A, = B, + B, A; + A,=B,+ B,, we find that 
every solution of (1) satisfies the system 
(A; + As) + (A; + Ay) + (Ap + As) + (Ag + Ay) 
= (B, + Bs) + (B, so B,) T (B, + Bs) +} (B, + B,) 
(A; + As) (Ay + Ag) (Ap + Ag) (Ap + Ay) 
= (B, + Bs) (B, + B,) (B, + Bs) (B, + B,). 


Therefore, from (7), we get the following identities 
(6ab + 6b”) + (a? + 8ab + 15b?) + (2a? + 10ab + 12b?) 
+ (3a? + 12ab + 21b?) 
= (2a? + 14ab + 20b*) + (3a? + 12ab + 9b*) 
+ (6ab + 18b?) + (a? + 4ab + 7b?) 
and 
(6ab + 6b?) (a? + 8ab + 15b*) (2a? + 10ab + 12b*) 
x (3a? + 12ab + 216?) 
= (2a? + 14ab + 20b%) (3a? + 12ab + 9b?) (6ab + 185?) 
x (a? + 4ab + 7b), 
by means of which we take particular solutions of the system (11). 
Thus, 
64+ 12+ 15+21=7+4+9+ 18+ 20 
6+ 21 = 124+ 15=7+20=9+ 18. 
6.12.15.21 = 7.9.18.20. 


IV 


The relation (8), for a= — a) + 4by, b =a, may be finally written 

Ag” + Agby — 3b”, 3 (ao” — agbp + bo"), ao” — 3agby — by?, 

X 3 (Gy? — agby + bo”) 
= 3a? — Saqho -+ bo?, do? + 3agby — bo”, 3ao2 — Apby — by”, 
ay” =e 5aybo oe 3b,y”, (n = Li 3. 5) 
whilst (7), for a= a)+ bo, b= ay) — bo, gives 

3ay” eins Sabo + bo’, 3 (3a9" ae 3ayby + bo”), 3a," supe Aybo ie bo, 
3 (39? — 3aqby + bo*) 
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= 9a," at Tagho + by’, 3a,” —_ Tagby -L 3b,", 9a,” — 1 lagb, 
a 3b,", 3a,” on Apdo —- bo” (n = 1, = 5). 


Adding the foregoing two relations, we take, after the omission of the 
index 0, the new relation 


a® — 3ab — b®, a? + ab — 3b*, 3(a?—ab+ 5%), 3(a@?—-a Ff 
+ b%), 3 (3a? — 3ab + b), 3 (3a? — 3ab + 5?) ; 
~ a? + 3ab — b®, a® — Sab + 3b?, 3a® — Tab + 3b?, 
3a? + ab — b®, 9a? — 1lab + 3b%, 9a? — Tab +B 
(n=1, 3, 5), (12) 


by means of which we get particular solutions of the system 





Ay, As As, As Ay As By Be, Bs By Bs Be (27 =1,3,5. 


Thus, for a=2, b= —1, we have 


nea ereinspeneinl im Ek = 


9, (—1), 21, 21, 57, 57 —3, 17, 29,9, 51,61, (@=1,3,5 ff 





that is 





3, 21, 21, 57, 57 = 1, 17, 29, 51, 61. (n = 1, 3, 5) 
V 
By interchanging the role of a and b we have, from (12), 

— a* + 3ab + b*, 3a* — Sab + b*, 3a* — Tab + 3b?, 

— a® + ab + 3b*, 3a? — llab + 9b?, a? — Tab + 90° 
= — a® — 3ab + b®?, —3a?+ ab+ b?, 3 (a? — ab + b?), 

3 (a — ab + b*), 3 (a? — 3ab + 36%), 
3 (a? — 3ab + 3b") @= i, 3,3 





Adding the above relation to (12) we take : 
3 (3a® — 3ab + 6), 3 (3a* — 3ab + b?), a? + ab — 36%, 

a? — Tab + 9b?, 3a® — llab + 9b, — a? + ab +364, 

3a? — Sab + b? ( 
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= 3 (a? — 3ab + 3b’), 3 (a? — 3ab + 3b’), —3a*+ab-+-b? 
9a* — Tab + b*, 9a? — llab + 3b?, 3a? + ab — b?, 
a® — Sab + 3b? (n= 1, 3, 5). 


This affords a solution of the system 








, Bi, Pe, Bag, Bays 320: , As = B, B, B,, B,,........ , Bs (n = 1, 3,5), 
depending upon two parameters. 
VI 
) } In order to obtain solutions of the system 
A,, As, As = B, Bo, Bs, By (n=1,3,5) (13) 
) we may proceed as follows: 
Starting with the system 
: M,, Me, Ms, My = Ny, No, Ns, Ny (14) 
and putting 
4(M, + M, + M;+ M,)=H 
» we may write the relation 
M, —H, Ms—H, M,;—H, M,—H 
= N, —H, N,—H, N,—H, N,—H, 
(m= 1, 2, 3, 5), ~— (15) 
: b where the sum of the first powers of the terms of each member is equal to 
4 zero. 
i Now, if 
5) ff H=N,, 
| we get 


Nz, — Ny, Ms — Ne, Ma — Nz 


— Ng — No, N, al Ng, Np ‘a M,, 
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If one applies this method in the particular relation 
2, 45, 76, 85= 1, 52, 65, 90, 
in which 
4(2 + 45 + 76 + 85) = 52 = Ng, 
one finds that a solution of the system (13) is afforded by 


52 — 1, 76 — 52, 85 — 52 " 65 — 52, 90 — 52, 52 — 45, 


52 — 2, (n = 1, 3, 5) 
that is 
51, 24, 33 = 13, 38, 7, 50, (nm = 1, 3, 5) 
where 
51 = 13 + 38. 


Then we shall confine our attention to the system (14), under the 
restriction 


N. = $(M, + Mz + Mz + My). 
For this, we consider the system 
2v, z+x—y, 2z—2(v +x), z+x+y 
=v, z z+, 2z — (v + x), 
in order to find a solution of (14), in which the supplementary conditions 
M, = 2N;, M. + M, = 2Nz3 
exist, as in the case of the above example. 
From the system under consideration we get the equivalent system 
3v2 + 3vux + 2x? -—z(2v+x+y?=0 
zx (7x — 3z + 12v) + 2zv (3v — 2z) — Tox (v + x) 
— 2x + 2(z+ x) y?=0, 
from which, eliminating y*, we find that 
22 + 6x? + 13v? — 5zx — 10zv + 13xv = 0, 
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Starting from the particular solution z= 2, x= —1, v =1 of the 
last equation, we put z= 2r+u, x=-—r+w, v=r._ Then the above 
equation gives r as a rational function of u, w, so that 

z = 3u® — uw + 12w? 
x = — u? + b6uw + 3w? 
v = u? — Suw + 6w? 
and the first equation of the equivalent system becomes 
y? = 4 + 8u8w — 1lu?w? — 18uv?. (16) 
Therefore (14) may be written 
2u? — 10uw + 12w?, (2u? + Suw + 15w?) — y, 
6u? — 4uw + 6w?, (2u? + Suw + 15w*) + yp 
i = uy? — Suw + 6w*, 3u2 — uw + 12w?, 2u® + Suw + 15w? 
6u? — 3uw + 15w? 
where the integers u, w, y satisfy the Eq. (16). 

Methods for resolving such equations as (16) have been given by several 
writers. 

Hence, we conclude as follows: 


< Particular solutions of the system (13) are given by means of the 
4 equalities 
2u? + 4uw + 6w?, 3u? — 3uw — 6w?, (— u? + 6uw + 3w?) + y 


= — y2+ 6uw + 3w?, 3u? — 2uw + 3w?, u? + 9uw, 
(w2 — 6uw — 3w?) + y, (n = 1, 3, 5) 





where the integers u, w, y satisfy the equation (16), that is to say 
y? = u(u — 2w) (u+ w) (u + 9w). 


Also, from (15) and under the above restriction, we get 


. — (u? + 9uw), — u? + 6uw + 3w? — y, 3u? — 3uw — 6w?, 
: — u* + 6uw + 3w? + y 

= — (2u? + 4uw + 6w?), — u2 + 6uw + 3w?, 

342 — 2uw + 3w?, 
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where 

mu=l, 2 3, 53> 
Thus, for 

u=2, w=—3, y= 20, 

we get 
38, — 24, 7= — 13, 51, — 50, 33, | 

or 
MM, 33, S1=7, 13, 38, D (mn = 1, 3, 5) 

and 
50, (— 33), (— 24), 7 = — 38, (—- 13), 51 (m = 1, 2, 3, 5) 


For u= 81, w= 40, y = 2079 we get finally 


5893, 6001, 11907 121, 5200, 6586, 11894, 
(a= 1, 3, 5) ee. 


Here we remark that by means of the relations 








m 





R,, Re, Rg = Uy, Un, Us, Uy (m = 1, 2, 3, 5) 
—R,=R,+R,, 
we get 
(R, + 0)? + (R, + R.)? + (R, + R;)? 
= (U, + Us) + (Uy + Us)? + (W, + U,)? 
(R, + 0) (Ri + Re) (Ri + Rs) = (Wj +U) (W4+Us) (U4 +U,). ; 
Thus, from : 
51, (—13), (—38) = 50, (— 33), (—24), 7, 
we have i 


512 + 38 + 13? = 172 + 26% + 572 
51.38.13 = 17.26.57. 


And, to finish, we note that a two-parameter solution of (13) is afforded 
by the values 
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A, = @ + 2a‘b — 2a*b? — Sa*b* — 2ab4 
A, = ath + 3a°b? — 3ab4 — b5 
A, = a + atb — a®b? + a*b® + 3ab* + BF 
B, = a + ath — 2a*b? — 4a*b* — 4ab* — bb 
B, = 2a*b + 3a*b? — a*b* — abt 
B, = —a*b + 5a*b* + 4ab* + b 
B, = a + 2ath — a®b? — 4a*b? — ab 
(where A; = B,; + B,) a and b being arbitrary integers. 
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(Department of Chemistry, Osmania University, Hyderabad, A.P.) 


Received July 4, 1958 


THE condensation of 4-chloro-o-phenylenediamine with benzaldehyde and 
p-nitrobenzaldehyde under different conditions was studied by Fischer and 
Limmer' who reported the formation of dianils, which could be converted 
into the corresponding 1: 2-disubstituted benziminazoles. The position of 
the chloro group (5 or 6) in the benziminazoles was not ascertained by these 
workers. In a recent publication,” the present authors established by un- 
ambiguous synthetic methods that the chloro group occupied position 5. 


A systematic study of the condensation between 4-chloro-o-phenylene- 
diamine and two moles of seventeen aromatic aldehydes (other than ben- 
zaldehyde*) containing alkyl, hydroxy, alkoxy, halo, nitro and amino sub- 
stituents has been made in the present investigation with a view to finding 
out the influence of substituents in the aldehydes as well as of the chloro 
group of the diamine on the nature and relative yields of the reaction pro- 
ducts. The condensations have been carried out following the modified 
Hinsberg’s procedure,* and the results obtained are summarised in Table I. 


The reaction has been found to yield mono (I) and di (II) substituted 
benziminazoles as in the other condensations.** On account of the pre- 
sence of low-melting resinous materials formed in the condensations, much 
difficulty has been encountered in processing the reaction products. 


Cl. J\_ is Be Ky, 
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CH,—R H 
Ir I 
It was shown previously in the case of condensations with o-phenylene- 
diamine* and 4-methyl-o-phenylenediamine,* that mono and di anils are 
possible intermediates in the formation of mono and disubstituted benzi- 
minazoles. It was also shown that electron-releasing groups like hydroxy 
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TABLE I 


i Products of condensation of 4-chloro-o-phenylenediamine with R.CHO 














No. R- % yield of (II) % yield of (1) 

1. p-tolyl = x aa eta 2 
2. o-hydroxyphenyl Bh i a 66 25 

3. m-hydroxyphenyl a a Sa 95 

4, p-hydroxyphenyl ie a = 95 

H 5. m-methoxy-p-hydroxyphenyl - si 95 
, 6. 3:4-methylenedioxypheny]. . ie =: 73 13 
7. p-methoxyphenyl vy ~ ‘a 78 15 
1 8. o-chlorophenyl ics oe - ies 45 
9. p-chlorophenyl - ss ” 30 4 
10. 2: 4-dichlorophenyl - ey “se “ 42 

i 11. p-bromopheny! A "3 : 31 
i 12. o-nitropheny] a a ey ne 50 

13. m-nitrophenyl - si “a 87 
14. p-nitrophenyl oe - oi 68* 13 


15. 2-hydroxy-5-nitrophenyl} .. 
16. p-dimethylaminophenyl .. - ¥s 87 
17. 1-naphthyl ei os on 88 





* Compound reported in literature. + Dianil isolated in 95% yield. 


and amino present in the aldehydes facilitate the exclusive formation of 
disubstituted benziminazoles, while electron-withdrawing groups like halo 
le and nitro favour the formation of considerable amounts of mono deriva- 
| tives. The results presented in this paper further corroborate these views. 
i The isolation of a stable dianil in the reaction with 5-nitrosalicylaldehyde 
: which could be converted to the corresponding benziminazole by refluxing 
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in acetic acid, lends support to the intermediary nature of dianils in the 
formation of disubstituted benziminazoles. 


The influence of substituents in the aldehydes on the course of the 
reaction seems to be affected to a certain extent by the chloro group of the 
diamine. The condensation of the diamine with salicylaldehyde yielded, 
in addition to the 1: 2-disubstituted benziminazole, a considerable amount 
of the mono derivative. o-Chloro- and o-nitro- benzaldehydes afforded 
exclusively monosubstituted benziminazoles in their reaction with this 
diamine whereas with other diamines** these aldehydes gave appreciable 
amounts of disubstituted benziminazoles also. These results indicate that 
the chloro group present in the diamine favours the formation of mono- 
substituted benziminazoles. p-Tolualdehyde gave a very small amount of 
the mono derivative and large amounts of high-melting products. Such 
high-melting and indefinite products have also been met with in the reactions 
of the diamine with halo and alkoxy aldehydes, on account of which the 
yields of benziminazoles have been rather low. 


EXPERIMENTAL 


All m.p.’s are uncorrected. The micro-analyses were carried out by 
one of the authors (C. V. R.). Unless otherwise mentioned, ethyl alcohol 
is the solvent used for recrystallisation. 


General procedure for the condensations——4-Chloro-o-phenylenediamine 
(one mole) and aromatic aldehyde (two moles) were dissolved separately 
in a slight excess of glacial acetic acid, and mixed together in a conical flask, 
when evolution of heat was noticed. The reaction was allowed to proceed 
for one hour with occasional stirring. If any crystalline solid separated 
by this time, it was filtered and the clear filtrate was added to excess of crushed 
ice with vigorous stirring. The compound that separated was filtered, and 
the filtrate was made ammoniacal, when a small amount of another com- 
pound was obtained. The compounds thus isolated were purified and 
characterised. 


2-(p-Methylphenyl)-5-(or 6)-chloro benziminazole was obtained by con- 
densation with p-tolualdehyde, from the ammoniacal solution. It crystal- 
lised in prismatic rods, m.p. 228° (Found: C, 68-8; H, 5-0; N, 12-1; 
CisHuNoCl requires C, 69-3; H, 4:5; N, 11-5%). 


1-(0-H ydroxybenzyl)-2-(0-hydroxyphenyl)-5-chloro —_ benziminazole — was 


obtained from the dilute acid solution of the salicylaldehyde condensation. 
The compound came out as needles from alcohol-acetone, m.p. 187° (Found: 
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C, 68-1; H, 4:8; N, 8-0; CopH,;N,0O.Cl requires C, 68-5; H, 4-3; N, 
8-0%). 


2-(0-Hydroxyphenyl)-5-(or 6)-chloro benziminazole was isolated from 
the ammoniacal solution of the above condensation. On purification, it 
came out as tiny rods, m.p. 273° (Found: C, 63-5; H, 3-9; N, 11-8; 
C,s3HgN,OCI requires C, 63-8; H, 3-7; N, 11-°5%). 


1-(m-Hydroxybenzyl)-2-(m-hydroxyphenyl)-5-chloro benziminazole was 
obtained by condensation with m-hydroxybenzaldehyde mostly from the 
dilute acid solution and in small quantity from the ammoniacal solution. 
The product on recrystallisation gave tiny rectangular rods, m.p. 247° (Found: 
C, 68-0; H, 4:7; N, 8:4; C.pH,;N,0,Cl requires C, 68-5; H, 4°3; N, 
80%). 


1-(p-Hydroxybenzyl)-2-(p-hydroxyphenyl)-5-chloro benziminazole was iso- 
lated from p-hydroxybenzaldehyde condensation, as in the above case. It 
crystallised in rectangular rods, m.p. 230° (Found: C, 68-3; H, 4-0; N, 
8:0; CpH,;N2O.Cl requires C, 68-5; H, 4-3; N, 8-0%). 


1-(m-Methoxy-p-hydroxybenzyl)-2-(m-methoxy-p-hydroxyphenyl)-5-chloro- 
benziminazole was prepared by condensation with vanillin and isolated as 
in the above two condensations. It was obtained as needles, m.p. 213° 
(Found: C, 64:0; H, 5-2; N, 7:0; C..H,ygN,0,Cl requires C, 64-3; H, 
4-6; N, 6°8%). 


1-(3 : 4-Methylenedioxybenzyl)-2-(3 : 4-methylenedioxyphenyl)-5-chloro ben- 
ziminazole-——In the condensation with piperonal, a sticky solid separated 
from the dilute acid solution, which on crushing with a little methanol, 
yielded a granular solid. The solid crystallised in rectangular rods, m.p 
171° (Found: C, 65:6; H, 4:2; N, 7-0; C..H,;N,O,Cl requires C, 65-0; 
H, 3-7; N, 6°9%). 


2-(3 : 4-Methylenedioxyphenyl)-5-(or 6)-chloro benziminazole was isolated 
from the ammoniacal solution of the above condensation and was obtained 
as leaflets, m.p. 205° (Found: C, 61-5; H, 3-8; N, 9°8; C,,H,N.O,Cl 
requires C, 61-7; H, 3:3; N, 10-3%). 


1-(p-Methoxybenzyl)-2-(p-methoxyphenyl)-5-chloro _benziminazole_ sepa- 
rated from the dilute acid solution of anisaldehyde condensation and was 
purified by recrystallisation, first from alcohol and then from petroleum 
ether, prismatic rods, m.p. 125° (Found: C, 69-2; H, 5-2; N, 8-0; 
Cx. HyyN,O.Cl requires C, 69-7; H, 5-0; N, 7°4%). 
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2-(p-Methoxyphenyl)-5-(or 6)-chloro benziminazole was obtained from 
the ammoniacal solution of the above condensation and on purification gave 
plates, m.p. 180° (Found: C, 64:7; H, 4:6; N, 10°9; C,,H,,N.OCI 
requires C, 65:0; H, 4-3; N, 10°8%). 


2-(0-Chlorophenyl)-5-(or 6)-chloro benziminazole——A sticky mass was 
obtained from the dilute acid solution of the condensation with o-chloro 
benzaldehyde, which on crushing with petroleum ether gave a granular solid. 
The solid crystallised in rectangular rods, m.p. 125° (Found: C, 59-0; H, 
3-5; N, 10°6; CysHgNeCl, requires C, 59-3; H, 3-0; N, 10°6%). 


1-(p-chlorobenzy!)-2-(p-chlorophenyl)-5-chloro benziminazole was obtained 
from the dilute acid solution of p-chlorobenzaldehyde condensation. The 
crude compound was purified by recrystallisation from alcohol and then 
from petroleum ether yielding needles, m.p. 150° (Found: C, 62°3; H, 
3-7; N, 7:5; CspHisNeCl, requires C, 61-9; H, 3-4; N, 7-2%). 


2-(p-chlorophenyl)-5-(or 6)-chloro benziminazole was isolated from the 
ammoniacal solution of the above condensation and on recrystallisation gave 
clusters of prismatic rods, m.p. 225° (Found: C, 58-8; H, 3-4; N, 10-5; 
CisHsN2Cl. requires C, 59-3; H, 3-0; N, 10-6%). 


2-(2: 4-Dichlorophenyl)-5-(or 6)-chloro benziminazole separated from the 
glacial acetic acid solution in the condensation with 2: 4-dichlorobenzal- 
dehyde. The compound crystallised in tiny needles, m.p. 184° (Found: 
C, 52-1; H, 2:8; N, 9-5; C.3;H,N2Cl, requires C, 52-4; H, 2:4; N, 9-4%). 


1-(p-Bromobenzyl)-2-(p-bromopheny!)-5-chloro benziminazole——In the p- 
bromobenzaldehyde condensation, a solid changing to a black resin separated 
out from the dilute acid solution. It was repeatedly crystallised from alcohol, 
and finally from petroleum ether yielding the pure benziminazole, prismatic 
rods, m.p. 155° (Found: C, 50:7; H, 3-2; N, 6°4; C.oH,;N,CIBr, requires 
C, 50:4; H, 2:7; N, 5-9%). 


2-(o-Nitrophenyl)-5-(or 6)-chloro benziminazole was obtained by the 
condensation with o-nitrobenzaldehyde, from the dilute acid solution. The 
compound was washed with a little benzene and recrystallised, shining yellow 
plates, m.p. 187° (Found: C, 56-5; H, 3-5; N, 15-5; C,sHsN3O0.Cl requires 
C, 57:0; H, 2-9; N, 15-4%). 

1-(m-Nitrobenzyl)-2-(m-nitrophenyl)-5S-chloro benziminazole was isolated 
from the dilute acid solution in m-nitrobenzaldehyde condensation. It 
crystallised from benzene containing a small quantity of alcohol in 
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pale yellow needles, m.p. 193° (Found: C, 59-0; H, 3-5; N, 13-4; 
CoHisNsOuCl requires C, 58-8; H, 3:2; N, 13-7%). 


2-(p-Nitrophenyl)-5-(or 6)-chloro benziminazole separated from the glacial 
acetic acid solution in the condensation with p-nitrobenzaldehyde. On 
purification, the compound came out as clusters of yellow needles, m.p. 
257° (Found: C, 57-5; H, 3-3; N, 15:6; C,3HsN3;0,Cl requires C, 57-0; 
H, 2:9; N, 15-4%). 


1-(p-Nitrobenzyl)-2-(p-nitrophenyl)-5-chloro benziminazole was _ isolated 
from the dilute acid solution of the above condensation and on crystal- 
lisation yielded yellow prismatic rods, m.p. 227° (Fischer and Limmer! report 
m.p. 235°) (Found: C, 59-2; H, 3-7; N, 14-1; CzypHigN,O,Cl requires 
C, 58-8; H, 3:2; N, 13-7%). 


N: N'-Di-(2-hydroxy-5-nitrobenzylidene)-4-chloro-o-phenylenediamine sepa- 
rated out from the glacial acetic acid solution in 5-nitrosalicylaldehyde con- 
densation and crystallised from nitrobenzene in orange yellow needles, 
m.p. 260° (decomp.) (Found: C, 55-0; H, 3-2; N, 12-5; C.pH,3N,O,Cl 
requires C, 54-5; H, 3-0; N, 12°7%). 


1-(2-Hydroxy-S-nitrobenzyl)-2-(2-hydroxy-5-nitrophenyl)-5-chloro benzimi- 
nazole, prepared by refluxing the above dianil in excess of acetic acid for 
twelve hours, crystallised from benzene in shining yellow prismatic rods, 
m.p. 314° (decomp.) (Found: C, 55-0; H, 3-5; N, 13+2; C.9Hi3N,0,Cl 
requires C, 54-5; H, 3-0; N, 12-7%). 


1-(p-Dimethylaminobenzyl)-2-(p-dimethylaminophenyl)-5-chloro benzimina- 
zole—The reaction between the diamine and p-dimethylaminobenzal- 
dehyde was conducted at 100° as the reaction did not proceed satisfactorily 
at the laboratory tempreature. The benziminazole was isolated partly 
from dilute acid solution and mostly from the ammoniacal solution. It 
crystallised from benzene-acetone in shining pale yellow prismatic rods, 
m.p. 239° (Found: C, 71-1; H, 5-8; N, 14-5; C.4H2sN,Cl requires C, 71-2; 
H, 6:2; N, 13°9%). 


1-(1-Naphthyl methyl)-2-(1-naphthyl)-5-chloro benziminazole was obtained 
by the condensation with 1l-naphthaldehyde, from the dilute acid solution. 
The compound was washed with a little acetone and recrystallised, rectan- 
gular rods, m.p. 199° (Found: C, 80-5; H, 5-0; N, 6-8; C.sHigN.Cl 
requires C, 80-3; H, 4°5; N, 6°7%). 


A2 
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SUMMARY 


The condensation of 4-chloro-o-phenylenediamine with two moles of 
seventeen aromatic aldehydes under acidic conditions has been found to 
yield mono and disubstituted benziminazoles in varying proportions, 
The influence of substituents in the aldehydes on the mode of reaction and 
cyclisation has been found to be similar to that observed in the condensa- 
tions with o-phenylenediamine. 
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SELF-RECIPROCAL FUNCTIONS IN THE 
FORM OF SERIES 


By V. V. L. N. Rao 
(Banaras Hindu University, Banaras-5) 
Received June 24, 1958 
(Communicated by Dr. T. Venkatarayudu, F.A.sc.) 


1, THE object of this note is to find out a few Self-Reciprocal functions in 
the form of series. Following the usual notation, I denote a function f(x) 
as R,, if it is Self-Reciprocal for Hankel Transforms of order » so that it 
is given by. 


fx) = J In Oy) FO) Vxy ay, (1.1) 
where J, (x) is a Bessel function of order p. If u = — 4, f(x) is denoted 


as R, while for » = 4, f(x) is written as Rg. 


2. Dr. Willis? has obtained expansions for certain Integrals, in the 
form of infinite series. Further Dr. Brij Mohan! has obtained a class of 
kernels for Self-Reciprocal functions of Hankel Transforms. On apply- 
ing these kernels in the Integrals given by Dr. Willis, we find that the ex- 
pansions in series become Self-Reciprocal functions. Accordingly we show 
in the following lines that certain Self-Reciprocal functions can be repre- 
sented as infinite series. 


3. Dr. Willis? has shown that 


f fe e-Mitdy = 1 { f+ f° ® .. £29) +. soni 


m 
(3.1) 
Let f(x) be Re or Rs. Then according to Dr. Brij Mohan! the kernel 


er, (3.2) 
transforms 

R-(R;) into Rs (Re) 
Hence we find that 

f em f(x) dx, | (3.3) 


is Rs(R,) according as f(x) is Re (Rs). 
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Hence from (3.3) and (3.1), we conclude that 
e{fO + = a +259 + | (3.4) 


is Rs(R,) according as f(x) is Rg (Rs). 





4. Again, Dr. Willis? has shown that 


frorinvd = 3{70+ 250+ Pot... : 


(4.1) 
Let f(x) be Rs. Then by the definition of an R, function we find that 








F FQ) sin xy dy = f (@). (4-2) § 
Hence from (4.2) and (4.1), we conclude that 
1 @+ + 49+... (4.3) 


is R, if f(x) is Rs. 
5. Further, it has been shown by Dr. Willis? that 








JS fo boordy =F {tO- FF + EO 





2? /2x* : 
1.3.5 (0 : 
— Pt b+ oe} 6-1) § 


Let f(x) be R, or Rs. 





Dr. Brij Mohan! has shown that the kernel 


J 0 (x), (5 ° 2) 
transforms 


Re (Rs) into R; (Re). : 


Hence from (5.1) and (5.2), we conclude that 


f (0) 1- 3f (0) 1. 3. 5.f © (0) 
3{fO- a7ixt B7Ie — se tf 


(5.3) 


is R,(R,) according as f(x) is R, (Rs). 
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6. It has been further shown by Dr. Willis that 


Self-Reciprocal Functions in the Form of Series 





Jt hevdy=L0 + (EO 4 LO, 


1.3f@)_ 1.3.5) 
2 72x5 — 23 /73x 


Dr. Brij Mohan? has shown that the kernel 
xD Js uty (x); 

transforms R, into R,. 

Putting v= (4+ 1) we find that the kernel 


Seuss (x), 


transforms R, into R,,,1; and in particular we find that 


Ji (x), 
transforms Rs into Rg,». 


Therefore, we find that 


J f(y) Si @y) ay, 
is Rep if f(x) is Rs. 
Hence from (6.5) and (6.1), we further conclude that 


£@ LS FP O_ FO _ 1.3f/"O 
+ es tae Soe 


1.3 e+... 4. 





is Rep if f(x) is Rs. 
7. Again Dr. Willis? has also shown that 











§ 10) 22 ay =F 1 + {LOO _- LO 


(6.3) 


(6.4) 


(6.5) 


(6.6) 
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so that , 
Foy sinay 4, _1[ £20 _ f©) 
§ 10) 2 ay = 3 [5 r+ {LA O_ Lee 
f fo - 7 
+ vein }. (7.2) 
Dr. Brij Mohan! has shown that the kernel 
=, (7.3) is R 
transforms R. (R32) into Rg. (Re). 
Hence it follows from (7.2) and (7.3) that sR 
l[m [7° @ _{£[*@. [°@ 
<5 f) + | x3 + 5x0 — ed], (7.4) 
is Raya (Re) according as Re (Re,s). we f 
8. It has been further shown by Dr. Willis* that 
2 _va 1 f®@) ' whic 
J roxreo = YF { £0) + Fy 3 FO 
1 1.3 f®@ 
74 2 Po + ed whic’ 
2(fm ©, Z1\ f®() | 
+ x { ae —— Z3 x3 expa 
a4. f°? OO, : 
es es oere Gl) fF 
Putting x for x? and y for y? we find that " 
co @ 
LWY) sxyg, — 3 ee (0) : 
SEP ery 3104 79° Ee 
1 1 f%@) 7 
+2°$"o +... 
£90 , LIF) > isgi 
bie ‘i x *7 xt ~ ‘ , 
Z2 f® (0) 
+ yl + oo | (8.2) 
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Hence from (3.2) and (8.2), it follows that the series 


1 f‘» () 1i1f"@ 
ae {f+ 7: “7 + ae fo tet 
+4{F00 42 Zl +s f® O) 








x? is ¥ 
+ seoohs (8.3) 
is Re (Rs) if the function 
f(VyY) 8.4 
vy’ lie 
is R; (Rc). In particular if 
f(vy) = 1, (8.5) 
we find that (8.4) reduces to the well known function? 
1 
—- 8.6 
Vy (8.6) 
which is Rs (R,). Accordingly (8.3) also reduces to the function 
1 
Vx? (8-7) 


which, then becomes Rg (Rs). 


9. Finally it has also been shown by Dr. Willis* that the Asymptotic 
expansion of the function 





f eS, (xy) ay, (9.1) 
) 
is 
l y* } 13.3% 1 1*.3*.5* 1 
s* 7ix9* JE: pt 2 ties 
Putting »/2x for x, and y//2 for y, we find that the Asymptotic expansion of 
a J eTa on) dy, (9-3) 


0 


is given by 


2 a 
_ + vat == sl Z2 ” F578x3 v2 y Vi2y7 +. 


@. 4) 
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Also Dr. Brij Mohan? has shown that 


ex*/2 (9.5) 
is R,. Hence from (9.5) and (5.2) it follows that 


8 (x) = Jo f ev 359) dy, 9.6 


is Rg. 
Also, it has been shown by Willis® that 


co 


2 
f EY" Jy (xy) dy = 7 E*'* I, (5). 


LY) 


so that 
. 2 
f ey’/2 Jo (xy) dy = /2n ex’! 4, (7). (9.7) 
Hence from (9.7), (9.6) and (9.4), we conclude that the R, function 
-x2 2 
e+ Ip (3) (9.8) 
has its Asymptotic gin given by the series 
=| |* 1%,.3* 1 1*,3%.5% 1 
NE: + at FF Bat 73 et} 
(9.9) 


I offer my thanks to Dr. Brij Mohan for his help and guidance in the 
preparation of this paper. 
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LILAVATI-A NEW ANALOGUE COMPUTER FOR 
SOLVING LINEAR SIMULTANEOUS EQUATIONS 
AND RELATED PROBLEMS 


Part II. Design of Model II and its Application to the Solution of 
Secular Equations 


By E. V. KRISHNAMURTHY 
(Department of Physics, University of Madras, Madras 25) 


Received July 27, 1958 
(Communicated by Prof. G. N. Ramachandran, F.a.sc.) 


1. INTRODUCTION 


THE general principles involved in the design of a new analogue computer 
(Lilavati) for solving linear simultaneous equations were discussed in a 
previous paper (G. N. Ramachandran and E. V. Krishnamurthy 1958— 
referred to as Part I). It was shown there how, unlike other computers 
where loading effects introduce errors due to current drain and non- 
linearity, this computer could provide an exact analogy of the set of equa- 
tions to be solved, by making use of Ohm’s law for multiplication and null 
methods for comparison of voltages. Even though these were the methods 
used in Model I for setting up the analogy, the measurements of voltages 
and currents were only done by meters, which are in principle not capable 
of highest precision. If exact measurements are to be made, only poten- 
tiometric methods should be used and meters should be completely dispensed 
with. Quite apart from the question of accuracy, the potentiometric arrange- 
ment has also other advantages; 


(a) Considerable economy is possible in material and in number of 
operations ; 


(b) The circuit can easily be modified to solve secular equations; 


(c) Operations can be made both with a.c. and d.c. without changing 
the circuit appreciably. 


This paper primarily deals with the design of Model II of Lilavati, 
employing potentiometers and its applications. Still further improvements 
are possible if one makes use of different iterative methods for solving both 
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simultaneous and secular equations and if corresponding changes in cir- 
cuitry are made. A discussion of these is reserved for Part III. 


2. THE POTENTIOMETRIC ARRANGEMENT 


The essential idea, as was mentioned in Part I, is to measure both 
currents and voltages by means of a standard potentiometer, the latter 
directly and the former by incorporating a standard resistance in the current 
circuit and measuring the voltage developed across it. In building Model II, 
the main aim was to test the suitability of potentiometric circuits in the 
design of the computer and consequently all components were designed 
only to have an accuracy of two significant figures (or 1%). 


The potentiometric arrangement itself is shown in Fig. 1 (a) and it 
could be used for measuring any voltage between 0 and 1-1 Volts to an 
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Fics. 1 (a), 1(5) and 1(c). Potentiometric arrangement for measuring current and voltage. 


accuracy of 0:01 V. It consists of a two-decade resistance box (R), the 
current through which could be varied by means of a variable resistance 
R,. In order to standardise it, a standard Weston cell is connected across 
101 2, with a tap key K and a galvanometer G. If the key K is pressed 
and the rheostat R, is adjusted for no current in G, then the potentiometer 
is standardised. 
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By means of two 11-way switches A and B, whose 11 ways are connected 
respectively to the points a» to dy) and by to by of the two-decade resistance 
box, any voltage value between 0 and 1-1 V. can be obtained across A and 
B, according to the setting of the switches A and B. 


In order to measure the current, a standard resistance of 100 2 (r) is 
included in the current circuit and the potential drop across it is compared 
with that between the points A and B of the potentiometer as shown in 
Fig. 1(5). The switches A and B are adjusted for balance and their setting 
gives straightaway the currents in units of 0-1m Amp. The same method 
can be used to measure any voltage in the range 0 to 1-1 V. in units of 
0-01 V. [Fig. 1 (c)]. 

3. THe New CIRCUIT 


In order to solve a set of equations: 


QyX, + cs se ee 8 8 8 * th Oy Xp eB (1) 
Qo1X = . . . ° . . . ° . - aon Xn = b, 

(1) 
AnyX, a . . . . . . . . . - Ann Xn _ bn 


the circuit shown in Fig. 2 is used. 


Here the values of a,s are all fed in as resistances on 100 2 wire wound 
volume controls by using an auxiliary Post-Office box arrangement, exactly 
as in Model I [see §4 (ii), Part I*]. On the other hand, the voltages b, are 
made equal to the required value by tapping the appropriate voltage from 
the potentiometer by switches similar to A and B. This is done by having 
a set of switches A,, B,; Az, B,; - + * *An, Bn similar to the set A, 
B for obtaining each of the values b,, bp - + + + + by respectively. 
This is possible for, in the Gauss-Seidel iterative process only one equation 
is set up at a time. 


The arrangement is shown in Fig. 3. The various potential points 
Ay tO Ay and by to by, in the two-decade resistance box of the potentiometer 
are brought to separate 1 pole 11 way switches A, A, - + + + An and 


* Experience showed that use of coarse and fine adjustments, as in Model I, was unnecessary 
if only an accuracy of 1% was aimed at. Hence each coefficient was represented in a singie 
100 2 volume control. The P-potentiometers which vary the current x, were also 100 Q 
yolume controls, 
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B, B, .... Bn as shown in the figure. Thus the voltages b, to by can 
assume any value between — 1-10 V. to + 1-10V. in steps of 0-01 V. 


depending on the setting of the corresponding switches, which can be readily 
read off on dials. 


; 
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Fic. 3. The setting switches. 


The cross-connections needed for setting up one equation at a time 
are shown in Fig. 2. The connections there are shown as having been made 
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for the first row. By means of a selector, as described in Part I [§4, (iv),] 
such connections can be made successively corresponding to each row of 
the equation from 1 to nm. The currents x, to x, are adjusted, as before, 
for balance at each step. The value of the currents need not be read at 
intermediate stages, but they can be accurately measured by means of the 
switches A and B, when the iterative process has converged. 


Thus in Model II, a single potentiometer can not only be used to obtain 
the voltage settings corresponding to the values b, on the right hand side, 
but also for measuring the currents xs. Further, the technique of measure- 
ment is based on a null setting, so that it is capable of great accuracy in 
principle. 

In §2, the method of standardising the potentiometer was described. 
‘n principle, there is no need to standardise the potentiometer, if the same 
potentiometer is used for setting the voltages b, and for measuring the 
currents x;. This is so since the value of the current is obtained from a 
measurement of the voltage drop over a known resistance, so that both the 
voltage values (b,) and the current values (xs) are affected in the same pro- 
portion if the current through the potentiometer is different from the 
standard value. The equations (1) are still true, since the resistance values 
ars temain unchanged. Of course, this argument is valid only if the current 
through the potentiometer remains steady during the process of solving 
the equations. 


However, if a dry battery is used for feeding the potentiometer, the 
em.f. may gradually fall during operation, thereby altering the current 
through the potentiometer. In such a case, it is necessary to keep the 
voltage across the ends of the potentiometer steady, by periodically stand- 
ardising it against some reference voltage. 


If a.c. is used, the need for standardisation can be completely eliminated 
provided the sources of e.m.f. for all the circuits are derived from the same 
primary source by means of transformers (see §5, below). 


4. DETAILS OF THE MopeL II (Fic. 4) 
(i) Description of the model 


The components used for building Model II of the computer were 
almost similar to those of Model I. The 100 2 coefficient volume controls 
(a,;) are mounted on the left hand side of the top panel (to the left of the 
galvanometer) and the sign switches for these coefficients are located on 
the right hand side. The toggle switch, at the top left of the galvanometer, 
is the one that is used for removing the shunt. The switch to the right of 
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this is the off-on switch for the voltage potentiometer. At the left of the 
lower panel are the on-off switches for the three current circuits and the 
plus-minus switches for currents xs. Below these are seen the P-potentio- 
meters P,, P, and P,; (with white knobs), for varying these currents. 


The multipole setting switches for b, (A;, B,; A», Bz; Ag, Bs) are 
located on the right of the lower panel, along with their sign reversal switches, 


The two knobs at the bottom are the setting switches A, B for measur- 
ing the currents xs. 


The selector switch S is at the top middle of the lower panel. The 
galvanometer is brought to the Wheatstone’s network in the first position 
of this switch (marked W) and in the second position (marked S) the gal- 
vanometer is brought to the standardising circuit. In the next three posi- 
tions of this selector (marked 1, 2 and 3) the iterative steps are carried out. 


The two-decade resistance box is fixed to the left-side of the computer. 
The 2 pole 9 way selector for bringing each of the nine a;s to the fourth 
arm of the Wheatstone’s bridge and the variable resistance (R’) that con- 
trols the current through the voltage potentiometer are mounted on the 
right side of the computer (not seen in the photograph). 


(ii) Auxiliary circuits 


(a) Circuit for setting resistances—The auxiliary circuits used in 
Model II were similar to those in Part I. 


The same two-decade resistance box R which is used for the poten- 
tiometer forms the third arm of the Wheatstone’s network. A D.P.D.T. 
switch located on the right side of the instrument enables one to isolate this 
box from the Wheatstone’s network while it is being used as a potentiometer. 
When it is switched off, the two contacts K,, K, in Fig. 5(6) are broken. 
A 2 pole 9 way selector switch brings each one of the volume controls (a;s) 


to the fourth arm of the Wheatstone’s network. The galvanometer is 


brought to this network by turning selector S to the position W. 


(b) Standardisation of the potentiometer.—In the last section, it was 
mentioned that while using a dry battery, one should keep the current through 
the potentiometer constant. For this purpose, a standard Weston cell is 
connected across 101 2 of the box, through the galvanometer, which may 
be done by bringing the selector switch to the position marked S. The 
variable resistance R, at the right side is adjusted for null reading of the 
galvanometer. The standardisation can be checked whenever necessary. by 
bringing the selector S to this second position (S). 
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(c) The sign switches——The sign switches for affixing plus or minus 
signs to Xs, by and ays are exactly similar to those mentioned in Part I §4, (i). 


(d) The equation selector—The equation selector used is an 8 pole 8 
way ganged switch. The function of this switch will well be understood 
from Fig. 5(a). The terminals marked C, D, E, F, G, H, I and J are the 
8 poles of the switch and cj, di, e, fi, gi, Ai, i; and j; (i= 1 to 8) are the 
eight ways corresponding to each one of these poles. 


For example, if this selector is set at the 4th position, since each pair 
of the poles (C, D), (E, F) and (I, J) are shorted, the points (c, and d,), (e, 
and f,) and (i, and j,) will be connected respectively, and the galvanometer 
will be connected between the points g, and hy. 


It is easy to see from the figure that at each of the settings 3, 4 and 5 
it connects up a particular row and also connects the galvanometer for 
balancing against the corresponding 5,. The iterative process is thus 
carried out in the 3rd, 4th and 5th positions of this switch (marked 1, 2, 3 
in Fig. 4). 

Fig. 5(b6) and 5(c) show how the same galvanometer is brought to the 
Wheatstone’s network and to the standardising circuit, in the first and 
second positions respectively of the selector switch. 


The 6th, 7th and 8th positions of this selector are used for reading off 
the currents x,, x, and x, [Fig. 5 (d)]. 


(e) Current measurements.—A standard 100 2 resistor is included in 
each one of the circuits and the ends of each one of these resistors can be 
brought for comparison against points A and B of the voltage potentio- 
meter (R) successively, in the 6th, 7th and 8th positions of the selector 
[Fig. 5(d)]. The switches A and B could then be adjusted each time for 
balance in the galvanometer and the current readily read off on the dials 
A and B. 


(f) The galvanometer—A 0-500u#A ammeter with a resistance of 
100 2 is used as the galvanometer with a shunt of 100 2, which could be 
removed when needed by means of a switch. The sensitivity of this gal- 
vanometer was found to be sufficient to get an accuracy within 1%. 


5. OPERATION WITH A.C. 


As mentioned earlier, the arrangement described in this paper can be 
used with an alternating voltage supply which would replace all the batteries. 
Assuming that there is no inductive effect, there will not be any phase change 
introduced due to the current flow in the various circuits. If the source 
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of alternating voltage for the various circuits, including the potentiometric 
circuit, are derived from a number of secondaries of the same transformer, 












































then the currents would all be either in phase or exactly out of phase accord- . o 
ing to the way in which they are connected to the secondary terminals. pu 
Therefore, if a commutator is incorporated in this circuit and it is reversed, ¥e 
the effect with a.c. is exactly equivalent to reversing the sign of the current } e 
with d.c. Consequently, all the arrangements used in the d.c. model for f 
changing the sign of xs, bys; and dys can be adopted, without any further m 
modification, for the a.c. model also. zs 
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In testing out Model II, with a.c., a 0-100A meter with a germanium 
crystal rectifier was used as a null indicator and was found to serve the 
purpose very well. A step-down transformer with a number of secondary 
coils (220 V.-6:3 V.) was used as the source of current for the various cir- 
cuits. Consequently the fluctuations in the mains voltage equally affect 
all the ciruits, so that no errors are introduced. There is no need therefore, 
for standardisation. In practice, it was found that, even when heavy fluc- 
tuations of the order of 100 Volts in the input voltage (150-250 V.) were 
introduced artificially by using a Variac, the computer worked quite satis- 
factorily, viz., the balance point was not upset by such fluctuations. 


Even though the wire wound potentiometers used were not strictly 
non-inductive, the error introduced was less than 1%. If one desires to 
have higher accuracy, it would be necessary to use non-inductive resistances, 
ie., a good resistance box, for each coefficient a;s. This would require 
n® resistance boxes, which would make the unit very costly (although it 
eliminates the need for a Wheatstone’s network). However, since only 
one row is connected up at a time, it is possible to devise an arrangement 
with only n resistance boxes, by having a series of n® setting switches, simi- 
lar to those used for setting the m voltages b,, from a single potentiometer. 
Such an arrangement is proposed to be built up and will be described in 
Part III. For this purpose an arrangement of resistances similar to the 
Kelvin-Varley potentiometer is required, rather than the ordinary two or 
three-decade resistance boxes. 


6. MODIFICATION TO SOLVE SECULAR EQUATIONS 


Of frequent occurrence in many fields is the problem of determining 
those values of a constant A for which non-trivial solutions exist to the 
homogeneous set of equations 


Qy3X1 a e e e e e e e . e - aAnXn = Ax, 
QyxX, + ee ee Axe 
‘ (2) 
QnyX1 + ee AXn 
or, in matrix notation 
AX = AX (3) 


Such a problem is known as a characteristic value problem; values of X for 
which nontrivial solutions exist are called characteristic values (also eigen- 
values or latent roots) of the problem. These values of A are also referred 
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to as the eigenvalues of the matrix A, i.e., the matrix composed of the co- 
efficients a,s. The corresponding vector solutions X (having components 
Xs) are known as the characteristic vectors (or eigenvectors) of the problem 
or of the matrix A. 


Equation (3) can also be written in the form 
(A — Al) X =0 (3 a) 


where I is the unit matrix of order n. This homogeneous problem p ossesses 
nontrivial solutions if and only if the determinant of the coefficient 
matrix, i.e., | A — AI |, vanishes or 
|A—~Al]=0 (4) 
Expanding this determinant, one finds that A is a root of an algebraic equa- 
tion of degree n, known as the secular equation. The n solutions Aj, d,, 
. An which need not all be distinct, are the characteristic values or latent 
roots of the matrix A. 
It is to be noted that if any one of these values of A are substituted in 
(2) and solved for x,...... Xn only (n — 1) equations will be independent 
so that only (nm — 1) ratios are relevant to the problem. 
A more general eigenvalue problem is to find the values of A which 
satisfy the matrix equation 
AX = ABX. (5) 
However this will not be discussed here. 


Model II can easily be modified to obtain the eigenvalues of equation 
(2) by using an iterative method described by Frost and Tamres (1947). But 
it cannot readily be adapted to solve the more general type of equations 
(5). The various methods of solving equation (5) and their relative merits 
as far as they can be mechanised are reserved for Part III. Here the aim 
is to show that a simple modification of the circuit discussed above can be 
used to solve equations of the type (2). 


We know that the left-hand side of equation (2) can be set up in analogy. 
If now this equation is rewritten in a slightly different form, viz., 


QQix—ADM+ ° 7 8 8 8 tt + Any, = 0 
Qnx, + ee 


AniX + ee ee ae > 2 + + + (an — A) xn = 0 
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it will be seen how a circuit analogy of the complete equation can be set up. 
It is only necessary to have a method whereby the value of the resistances 
representing the diagonal terms 4aj,,...... Ann can all be decreased (alge- 
braically) by the same magnitude. 


This could conceivably be arranged by using ganged potentiometers 
for these resistances: In that case, one method of solving (3) will be to 
start with arbitrary values x; = 1, X., x3 .... Xn for the various currents 
and A of A and then to make use of the following iterative process (Frost 
and Tamres, 1947). 


(i) Set up the first equation and change x, to x, for balance. 
(ii) Set up the second equation and change x, to x, “) for balance. 


(iii) Continue until the (7 — 1)-th equation is set up and adjust the 
n-th variable x, to x, ™ for balance. 


(iv) Set up the n-th equation and adjust A to A“ for balance. 


(v) Repeat the whole sequence until the quantities x2, x,...... Xn» 
A remain constant within the desired limits. 


However it is difficult to gang potentiometers accurately and so a dif- 
ferent method was adopted for Model II. In this method, an exact analogy 
of equation (2) is set up. Here the left-hand side is exactly the same as for 
the solution of simultaneous equations. On the right-hand side, however, 
we do not have the potentiometers giving voltages b,.... bn, but a single 
potentiometer (A) whose resistance can be set equal to A and the current 
through which can be made equal to any required value from x, to Xn. The 
latter adjustment can readily be made by incorporating a standard resist- 
ance r (100 2) in series with the A potentiometer and comparing the voltage 
developed across this with an equal resistance r (100 2) kept in series with 
each of the current circuits x, to x, (Fig. 6a) [the latter are already there 
in Model II, see §4, ii(e)]. However, it is not necessary to have a series 
resistance along with the A-potentiometer; for one can tap off the required 
resistance r (100 2) directly from the A-potentiometer as shown in Fig. 6 (0), 
and bring it for comparison, so that equalisation of the current can be made. 
The advantage of having this arrangement is that the setting switch A-B 
(which is used for measuring the currents if the potentiometer is standardised) 
can directly be made to serve this purpose, provided it is adjusted to read 
the value r (100 2). In order to set up the analogy of the right-hand side 
of all the equations (2), we need only a single potentiometer, through which 
the current flowing can be made different each time. So instead of having 
all the n(Ay — Bn) setting switches any one setting switch (say A, — B,) 
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may be chosen for the purpose, the current through which being made equal 
to any one required value by the above method. Also, the value of A can 
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6 (a) 6 (6) 
Fics. 6(a) and 6(6). Equalisation of current in A-potentiometer. 


be directly read off on the dials. These arrangements are easily made, if 
a slight modification is made in the connections to the selector switch §, 
shown in Fig. 5(a). The modified circuit is shown in Fig. 7. 


The terminals C, D, E, F, G, H, I and J are the 8 poles of the switch 
and cj, di, ei, fi, gi, M4, i; and jj (i= 1 to 8) are the 8 ways. We will now 
see how the currents are equalised using this selector. Since C and D, E 
and F, and I and J are shorted and the galvanometer is connected between 
G and H, if the switch is set to the third position the points (c, and ds), (e; 
and f;), (is and j,) are respectively connected and the galvanometer is intro- 
duced between g3; and fA. This serves to make the current through the 
A-potentiometer equal to x,, the current through the first circuit. Simi- 
larly, at settings 5 and 7, one can make the current through the A-potcu- 
tiometer equal to x, and x;. At these same positions 3, 5 and 7, it is also 
possible to measure the currents x,, X2, X3 if the potentiometer is standardised. 
Thus the following will be the sequence of operations that will have to be 
made in order to solve equations (2) by the iterative process mentioned 
above. 


(a) Selector is set in position |.—The coefficients ays are all fed in as 
resistances using the Post Office Box arrangement and the signs of ays are 
also set. 


(b) Selector is set in position 2.—The potentiometer is standardised 
by adjusting R, [as in §4, (ii)]. 
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(c) Selector is set in position 3.—The current through the first circuit 
js measured by means of the setting switch A-B and adjusted to some 
standard value or unity (say x). 


(d) The setting switch A-B is set to read 100 2 and with selector in 
position 3, the variable resistance R, is adjusted for balance in the galvano- 
meter. Thus the current x, is made to flow through the A potentiometer. 


(e) The other variables xg and x, are represented as arbitrary currents 


x. and x3 need not be measured and in any arbitrary value for A, say 
\ is set on the switch A,-B,. 


(f) Selector is set in position 4.—Potentiometer P, is adjusted and 
x, is varied to x," for balance. 


(g) Selector is set in position 5.—The variable resistance R, is adjusted 
for balance making x,") to flow in the A-potentiometer. 


(h) Selector is set in position 6.—Potentiometer P, is adjusted and 
x3 varied to x3") for balance. 


(i) Selector is set in position 7.—The variable resistance R, is adjusted 
for balance making x," to flow in the A-potentiometer. 


(j) Selector is set in position 8.—The setting switch A,-B, is adjusted 
for balance to get a better value for A, namely A”. 


(k) The selector is now brought back to position 3 and R, is adjusted 
for balance in galvanometer, thus starting a new cycle of operations. 


These cycles of operations are continued until the quantities x,, x, 
and A remain constant within the desired limits. Then the values of x, 
X3 are measured and A can be read off from the setting of the switch A,-B,. 


This iterative method presents no difficulties in converging to the largest 


and smallest eigenvalues. Once these eigenvalues are obtained, it is possible 
in principle, to force the convergence upon some other value by keeping 
the trial vector x orthogonal to both the largest and the smallest eigenvectors. 
But in practice, some difficulties in convergence to the intermediate eigen- 
values may be encountered unless a good approximation is available to 
start with. This difficulty is felt even in numerical methods involving 
this particular iterative process. The method very largely depends on the 
experience of the computer and so automatic methods cannot be adopted 


(Murray, 1948). Fuller details on other iterative methods will be discussed 
in Part III, 
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The computer (Model II) can be present for solving simultaneous equa- 
tions or for solving the secular equations by making use of a 6 pole 2 way 
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Fic. 7. Modified circuit for secular equations. 
switch T (Fig. 8) (not seen in the photograph). When this switch is in its 
first way, it connects the points (4, and A,), (j, and B,), (4g and Ag), (jg and 
B,), (hg and Ag), (jg and Bs) so that the computer can be used to solve simul- 
taneous equations. When the switch T is in the second way, the computer 
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is ready for solving secular equation, since all the points h,, hg and hg are 
connected to A,, and j,, jg and jg are connected to B,, as shown in Fig. 8. 


As a result of this modification, namely of the computer being able to 
solve either simultaneous equations or secular equations, the settings of 
the selector switch at which the iterative process for solving the former are 
made, are slightly different from those described in §4, ii(d). Instead of 
iterating in the settings 3, 4, 5 and measuring the currents xs at 6, 7 and 8, 
iteration is now performed at settings 4, 6 and 8 and currents are measured 
at 3, 5 and 7. 


The computer described here can handle only a third order secular 
equation. However there is no essential difficulty in constructing a larger 
computer based on these principles by making use of switches having larger 
capacity. 


7. SUMMARY 


This paper essentially deals with Model II of Lilavati (an analogue 
computer for solving linear simultaneous equations) where unlike in the 
first model, potentiometric methods are employed for measurements of 
current and voltage. The incorporation of the potentiometer has con- 
siderably simplified the number of operations required and has also helped 
in successfully modifying this computer to solve characteristic value prob- 
lems. The operational details of this computer are also described in this 
paper. 
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ABSTRACT 


Some polishing experiments have been carried out on copper anodes 
in a hyperbolic cell designed by Gilmont and Walton, using orthophos- 
phoric acid as the electrolyte. The results obtained have been compared 
to those obtained in similar experiments in a Hull cell. It has been 
found that very similar bands of different reflectivity and polishes are 
found to form in both the cells. These bands shift with time and a 
study of such displacements has been made. The results are briefly 
discussed. 


INTRODUCTION 


OF late various workers! have tried to design and work out the theory of 
different types of electrolytic cells which would give a uniformly varying 
current density along an electrode assuming that there is no polarization, 
a condition normally not attainable in practice. The use of the Hull? 
cell for studying electro-plating phenomena is well known. The construc- 
tion of the Hull cell is fairly simple, but a theoretical and mathematical 
analysis of the current distribution is not so easy. Although recently 
Skwirzynski and Huttly* have given a mathematical solution of the current 
distribution for a trapezoidal cell very much similar to Hull cell, yet they 
found that the current distribution is linear only on about ¥ of the electrode, 
there being a singularity at the high current density end. According to 
Kasper,’ the proper choice of electrode geometry should give a sufficiently 
linear current distribution with simple mathematical solutions which would 
eliminate the need for emperical methods. Recently, Gilmont and Walton,5 
on the basis of Kasper’s ideas, have designed a new cell, usually referred 
to as hyperbolic cell, in which the current density function along an electrode 
is almost linear and can be predicted by a simple equation. The correct 
geometrical shape of the new cell was obtained by applying the fundamental 
principles of the theory of electro-potentials, according to which the current 
density lines are orthogonal to the equi-potential surfaces, 
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Briefly stating, the hyperbolic cell (Fig. 1) as designed by Gilmont and 
Walton consists of 4 sides, two of which are straight and equal intersecting 
at 45°. The other two sides are rectangular hyperbolas orthogonal to each 
other and each in turn orthogonal to one of the straight sides. One of the 
straight sides and its opposite curved surface are used as electrodes. The 
straight side is used as the anode while the curved side as the cathode. 
Recently, we have reported*®? some of our experimental findings on the 
electrolytic polishing of copper in a Hull cell using orthophosphoric acid 
as electrolyte. There we observed highly different finishes on one and the 
same specimen, with sharp contrasts and discontinuities between these 
finishes or bands. Also, it was found that two almost equally polished 
bands were always separated by a small unpolished band. The position 
of the various bands was also found to shift with the time of polishing. As 
no definite and clear insight into these phenomena is so far available and 
as the current distribution in the Hull cell is rather complex, it was consi- 
dered worth while to compare these observations with similar experiments 
using the hyperbolic cell. The results of these comparative experiments 
in a hyperbolic cell are reported below. 
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In order to facilitate comparison, the hyperbolic cell was also made 
of perspex and the electrolyte used was of the same concentration (756 g/1 
of H;PO,) as employed in our Hull cell studies. The experimental details 
and copper (high purity 99-98%) electrodes were similar to those used in 
our previous study. 


RESULTS AND DISCUSSION 


Figs. 2, 3, 4 and 5 show the various effects produced on copper electro- 
polished in the hyperbolic cell for varying times. It will be seen that in 
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all the figures there are 5 well-defined regions showing different reflectivity, 
colour, structure etc. The same are briefly described and shown in Fig. 6. 
The above results were found to be quite reproducible both as regards their 
nature and their relative position. The bands are also found to be slightly 
curved at the top towards the high current density end and straight at the 
bottom. The most interesting observation was that although this type of 
cell is known theoretically to have a very linear current distribution, yet, 
the two well polished regions (II and IV, Fig. 6) were always found to be 
separated by a small unpolished or frittered region (III), exactly as was the 
case in Hull cell experiments. 
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The shift in the position of bands with varying time also similarly com- 
pares with that observed with Hull cell, i.e., as the polishing time increases, 
the different bands are found to shift towards the low current density end 
upto a certain extent of time after which they tend to retrace and shift back 
towards the high current density end. For example, during polishing at 
a total cell current of 6 A, the bands shift towards the low current density 
end upto a period of about 6 min. If the polishing time is further increased 
beyond 6 min. the shifts are in the opposite direction (see Figs. 4 and 5) as 
the time of polishing increases. Similarly, for Hull cell, it was found that 
at a total cell current of 4 A, the retracing of the band-shifts occurs after 
about 8 min. 


Thus, we see that on the whole the results obtained both in a hyperbolic 
cell and in a Hull cell resemble each other quite closely. The full logic of 
these observations is not yet clear, It appears that for understanding the 
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above phenomena, we must look for a detailed knowledge of (a) the pro- 
gressive effect of the anodic reaction on the composition of the film of liquid 
near the anode, (b) the formation of solid films on the anode and their elec- 
trical properties, (c) the effect of the resulting intense polarization on the 
current density pattern on the anode as a function of time, (d) the current 
density values necessary for getting polishing effect under the particular 
hydrodynamic conditions obtaining at the anode, and (e) short-circuiting 
effect by the anode on the regions of different potentials and the consequent 
effect of current density pattern. 


It appears that the actual conditions existing in such cells are very dif- 
ferent from those initially assumed for working out the cell design. In 
electropolishing, in particular, there is a high degree of polarization due 
to the production of poorly conducting anodic films. This leads to com- 
plete deviation of the current distribution on the electrode, from the theore- 
tical function. Since, moreover, the polarization is a function of time and 
depends on the hydrodynamic conditions, the whole pattern of current 
density distribution would considerably change with time. 


No detailed study of the above factors has so far been reported and 
further work is required for a correct elucidation of our observations. 
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1. INTRODUCTION 


THE first theoretical investigations of the diffraction of light by superposed 
supersonic waves in liquids and solids date from 1938. Theories were 
then developed by Bergmann and Fues,! Nagendra Nath? and Nagabhushana 
Rao* and were mainly based on Raman and Nath’s preliminary theory.‘ 
The considered sound waves have frequencies which are incommensurate 
and they do not present phase differences, so that the expressions for the 
intensities of the spectra derived by those authors only have a limited use. 


Experiments showing combination spectra were performed by Berg- 
mann? (in liquids), Hiedemann and Hoesch® (in solids), Bergmann and Fues! 
(in liquids and solids) and Govinda Rao’ (in liquids). None of those ex- 


perimental studies gives quantitative measurements of the intensities of the 
obtained spectra. 


In 1949 Ramachandra Rao® set up a theory of the diffraction of light by 
two superposed parallel ultrasonic waves, the wavelength of the second 
being half the wavelength of the first (or the ratio of the frequency of the 
first to the frequency of the second being 1:2; there also exists a definite 
phase relationship between the two sound waves. Using the method 
of Raman-Nath’s simplified theory and confining himself to a perturbation 
showing no time dependence, he comes to the interesting conclusion that 
the obtained spectra must be asymmetric with respect to the central image. 
This feature has been confirmed experimentally by the same author. 


Those calculations were extended by Murty® to the case of two sound 
waves, the frequencies of which are generally in the ratio 1:n. This theory 
which is still based on Raman-—Nath’s elementary theory leads to the fol- 
lowing conclusions: 


* Geassocieerde van het Nationaal Fonds voor Wetenschappelijk Onderzoek van Belgié. 
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When the ratio of the frequencies is 1: (27+ 1) the intensities of 
opposite orders are equal so that the pattern is symmetric; 


When the ratio of the frequencies is 1:2/ the intensities of opposite 
orders are different and the pattern is asymmetric. 


Recently the interest in the problem of the diffraction of light by 
superposed parallel supersonic waves was again revived by the study of the 
distortion of finite amplitude sound waves by Zankel and Hiedemann”: 
the distorted wave is resolved into its Fourier components, so that the latter 
can be considered as superposed waves. 


The purpose of the present paper is to expose the general theory of the 
diffraction of light by two superposed parallel supersonic waves, one being 
the n-th harmonic of the other. In the second section a system of differ- 
tnce—differential equations for the amplitudes of the different spectra is 
derived. In §3 the special case p = A?/uou,A* *? = 0 is considered and the 
system of difference-differential equations is solved exactly. In the follow- 
ing section it is shown that in the case p# 0 the diffraction pattern is always 
asymmetric with respect to the zero order, as well for even as for odd values 
of n. Finally (§§ 5-7) a series solution of the general system is given for 
n=2 and n= 3. In the latter case coefficients are calculated up to the 
term giving the asymmetry effect in the intensity. 


It is a great pleasure for us to pay tribute to Prof. Sir C. V. Raman 
for the fundamental theories for the diffraction of light by supersonic waves 
he has established and which are characterized by the simplicity of the basic 
ideas underlying them. Both his preliminary and generalized theories were 
the first to give a satisfactory explanation of the diffraction patterns contain- 
ing a great number of spectra and formed the basis for many theoretical and 
experimental studies, the present investigation included. 


2. ESTABLISHMENT OF A SYSTEM OF DIFFERENCE-DIFFERENTIAL EQUATIONS 
FOR THE AMPLITUDES OF THE DIFFRACTED LIGHT WAVES* 


Let us consider a parallel beam of light of frequency v and wavelength 
in vacuum A, passing through a liquid column, disturbed by two superposed 
parallel progressive supersonic waves, one of them being the n-th harmonic 
of the other. We further suppose that the directions of propagation of the 


* Independently and using the same method, Zankel" has obtained a system of difference- 
differential equations in the case where a complete set of harmonics, without phase differences, 
are contained in the sound wave. For this system, which is a generalization of our equation (9), 
he has found the exact solution for p = 0. If only one harmonic, containing a phase constant, 
is present Zankel has also obtained independently the exact solution (17) of Section 3, in the case 
p=0. 
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light and sound waves are perpendicular to each other. We shall put the 
z-axis along the direction of the incident light and the x-axis along the direc- 
tion of the sound waves. 


In order to obtain the differential system for the amplitudes of the 
different spectra we shall follow the method of Raman and Nath’s general- 
ized theory.!? 


Writing for the electric field of the light 


E (x, z, t) = e***¢ (x, z, t), (1) 
¢ (x, z,t) satisfies the equation 
Da\2 
4$=—({) #08; Q) 


p(x, t) is the refractive index of the disturbed liquid and is given 
by the expression 


a ae 2m (vt - is) 


+ pn sin 27 (nv*t = a + 4) (3) 
where 


\* = wavelength of the supersonic waves in the considered medium; 
v* = frequency of the supersonic waves; 
fo = refractive index of the undisturbed medium; 


#4, = maximum variation of the refractive index of the disturbed medium 
due to the fundamental tone of the supersonic wave; 


fn = maximum variation of the refractive index of the disturbed medium 
caused by the n-th harmonic of the supersonic wave; 


4 =phase of the second wave with respect to the first. 


According to the periodicity of ¢ (x, z, t) in x and ft, this function may 
be developed in a double Fourier series 


¢ (x, Z, t) = > z Srs (2) e2mir x/\* amisy *t (4) 
Taking into account that ¢(x + pA*, z, t) = (x, z, t— p/v*), p being any 
number, it follows that f;s(z)=0 when r#—s and f;s(z) #0 for 


r=- S. 
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Hence, 
¢ (x, z, t) = yr FA(2) e2nir C/A Q—2wiry*t (5) 


Substituting this series into the differential equation (2), equating the 
coefficients of equal exponentials on both sides, putting: 


fr= on tat mz }, (z) (6) 


into the result, and neglecting the second derivative of ¢,(z) with respect 
to the other terms, we finally obtain the following system of difference- 
differential equations, 


2 Br — oy rs — $ris) — On Grant? — drinti®) = iBr%p (7) 
with 
0 =, on = EM, B= ey, 8 = Ind (7a) 
or putting 
a... 2 B A? = (8) 





2a — (bra — brit) — On (Gr_ne-? — brine") = ipr?dy. (9) 








From (5) and (6) it is seen that the r-th order spectrum makes an angle 
6, defined by sin 6, = — r(A/A*), with respect to the z-axis, and undergoes 
a frequency change — rv*; its amplitude is proportional with ¢,, the factor 
of proportionality having modulus one. 


The boundary conditions are 
br (0) = dr9 (10) 


where 5,9 is the well-known Kronecker symbol. 


3. SPECIAL CASE: p=0O 


For great values of p+ and A* the right-hand side of equation (9) is 
negligible and the system becomes 


2a — (¢r-1 — Orsi) — On (¢r-ne“* one ¢rine”) = 0. (11) 
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In order to integrate the system (11) with the boundary conditions (10) 
we shall use the complex function method of Nagendra Nath-Wilson.? 
We shall put 


$f ( 
br) = a5; (FO) ar (12) 


where F() and f(y) are yet unknown functions of the complex variable 
n, and the integral is taken along a closed path in the complex 7-plane, 
encircling the origin once in the positive sense. According to the boundary 
conditions (10) it is required that the function F(y) must be 


F (n) = 1. (13) 


Substitution of the expression (12) into the equation (11) gives 


| an( nis __1 
Fo) = 3 (1-5) + Flare wag) (14) 
so that the solution becomes 


on [§ O—D]eo [°F (ne ab) 


Tr 
b 





br (= 55 ( () 
Xd. (15) 
Using the well-known Laurent series for the exponential functions!® 


exp.[ 5 (n— <) ]- yy) »or (16 a) 


p=— co 
exp. [48 (xe — 


(Jn (z) being the Bessel function of order n] 
in the right-hand side of (15) we obtain 


br (= 5, > bi Jp (6) Tq (ant) e198? (fy) gPtPa-T-1 dy, 


p=-CO g=-cOo 


=) | mt Tq (an€) n™4e*, (16 5) 


q=- co 


n 


The integral being equal to 27i5p rng, we finally have 





dr (€) = -£ ' Irma (§) Jq (ané) a 

















PS RAR AM me 
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which are in fact the solutions obtained by Murty, using the method of 
Raman and Nath’s preliminary theory, but written in a more condensed 
form. Calculating the intensities, Murty finds that in the case of odd values 
of n, the intensities of the spectra of orders + r and — r are equal, so that 
the spectra are symmetric with respect to the spectrum of order zero; in 
the case of even values of n the spectra are asymmetric with respect to the 
order zero, i.e., the intensities of the orders +r and —r are not equal 
(excepted for 5 = 7/2). 
4. PROPERTIES OF THE INTENSITIES IN THE GENERAL CASE (p+ 0) 


We shall now show that Murty’s conclusion for odd values of n will 
not hold in the general case. 


Let us suppose that the intensities of the orders r and — r should be 
equal, i.e., Ip—= Ly, then we should have either ¢,-=e'%¢., or 
by = eld *, 


(a) First possibility: 
or = eG ,, 


Introducing this substitution into the general equation (9) we obtain 
2 ae — [upg 8 (M-- 1) — py (YH - 1) 
— On [Ppyn © (Ye-n— Yr) E18 — h_p_n Ci (Yeon) C8] 
= (-2 oH + ipr*) 4 (18) 
On the other hand, changing r into — r in equation (9) gives 
2 oe ins — Kind ~ anion — dice 


= ipr*d_,y. (19) 
Equations (18) and (19) must be identic, hence, 
(1°) dy-/dé = 0 from which it follows that y; is a constant; 


(2°) ef (%--%r) = — 1 and ef(%u-%) = — 1, from which 
Yr =r, (20) 
(3°) ef (Yn) €-1 = — e's and ef (rH 71) CF = — E-18, 


from which taking into account (20), 


cos 8 = 0 or 8= (k + 4) 7 if nis even, 
Ad 
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sin 85=0 or 8= ka when n is odd. 
(b) Second possibility: 
be = err (8)b_p*. 


Making this substitution into the system (9) gives 


? db_,* a ‘( - *  ( —%7r) 
— [P¥ 941 8 (Yea te) — P¥_p_1 (Wr -Y)] 


— Gy [$*_rin e-7d e: (¥r-n— Yr) _ p*_,_net6 e (Yr+n—r)] 
ane +. dyr k 
- i(pr 2 E) $*_,. (21) 


Changing r into — r in the equation (9) and taking its complex conjugate, 
leads to 


y dp_,* * * * 16 * —i5 
dé — (P¥_p-1 — $* 741) — On (6*_7--ne”? — b¥_p:nO”) 


= — ipr*d_,*. (22) 


The identity of the second members of the equations (21) and (22) gives 
dy,|dé = r*p from which y, = pr?é + rm (when p= 0, then yp =rz). The 
first bracket in the left-hand side of equation (21) then reads ¢*_,_, e#?§(1+2") 
— $*_,,,e%8"-2"), Comparing this expression with the first bracket of the 
left-hand side of equation (12), it seems that their identity could never be 
realized. Hence, we may conclude in the general case, p#0: if n is even 


and 8#4(k+4)z or if n is odd and 8+ kz the intensities of the orders t 
and —r could never be equal. 


In the case p= 0, the same reasoning leads to the conclusion that the 
intensities of orders r and — r cannot be equal if m is even. When n is odd, 
the second possibility predicts that the intensities of orders r and — r could 
be equal and that the relation between ¢, and ¢_, should be ¢, = e!*"4*_, 


= (—1)" ¢,*. It is easy to verify that the solution (17) obeys this relation 
for odd values of n. 


Ramachandra Rao’s experimental observation’ of a symmetric spectrum 
for n = 3 could be explained by the fact that his experimental conditions 
are so that p<1 or that the phase angle 5 = kz. 
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5. SOME GENERAL CONSIDERATIONS ON THE SERIES SOLUTION OF 
THE GENERAL SYSTEM (9) 


co 
Introducing the series ¢,= ¥ A;~é* into the system (9), leads after 


k=0 
careful inspection to the conclusion that the series must have the following 


form: 
(a) m= 2] + 1, me = 0, 1, 2, ...... 
bo = 3 Age 
k=0 
— k 
P+ (21-1) m+1 == £1 DA + coliary mt1,k§ 
k=0 
+ k 
+ (2litymt2 = EM? DAs col mt2,k& 
k=0 
SP k 
$+ (2141) m+l = em © A + c2tvy m+ ké 
k=0 
lia Fp k 
b+(2liymtilay = Et DA + (2141) m+ (lin, kS 
k=0 
ly k 
$+ (ali m+i1+2) = §mt YS A+ (oli m+ le, k§ 
kzw0 
7 k 
b+ alan m+ cel) uo 2S A+ (alin m+ clin, ké 
k=0 
(6) 2 = 2], m= G@, I, 2, ...... 
~ k 
do == 2D Aoké 
k=0 
> k 
b+olm+1 = EM SA + otm41,k§ 
k=0 
PR 
P+olmte = §™? FAs a1mt2,kb" 
k=0 
mil k 
¢+olm+l = @m 2 Asoimtt,é 
k=0 
m+l k 
$+2lm+il) = EM" DY Asolmstil-y, kf 





k=0 


co 
== §™ ‘Z A +elmial, ke”. 
Lad 


t (23) 
} (24) 
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It may easily be verified that the lowest powers of é in those series are 
the same as those in the solution (17) of the system in the special case p = 0, 


6. ExpLicir CALCULATION OF THE INTENSITIES BY THE SERIES 
METHOD IN THE CASE n= 2 


According to the series expansions (24) we may write for n = 2, 


bo = Aoxe* 
k 


a, em 1 
d 9 = A 2m+i1,k 
t+2mt+1 Zm+ Im! +2m+1, gk 


k=0 
m-+1 eEm+1 


P+omte = = eH Get 1)! 





The factors a.™/2™*!m! and a,™+1/2™+1 (m+ 1)! are written for con- 
venience and are chosen as the corresponding lowest power coefficients of 
the series expansion of the solution (17) for p= 0. Substituting the series 
(25) into the equation (9) and equating the coefficients of equal powers of 


€ on both sides, we find the following relations between the coefficients 
Ark? 
(m+ k + 2) Atomsi, kus F Atom, ka + Ty “s i) 


a,” 


+ MA + om¥1, kuie*? = 4(m+ J) 4 +2mts,k- jet% 


=F Qm+1)?Asemtik  (m>1) 
(kK + 2) Asi, keu F Ao, kw + > Ate, k ss > Azi,KeF® 
x “Ass kiet® 


= F Ask (27) 


1 
(m +k + 2) Asames, ea F 5 Asomtikt4 Atom+s,k-1 


F (m+1) A+om, ky eF ® 


2 . 
* tm 3) A+om+ 4,k-et” 


=F (m+ 2)*Asemssk  (m>0) (28) 








re 


5) 


6) 
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2 
aa 


(k + 1) Ao ku a 4 A_,, ka + 4 Ay k-1 oe A_,, k-1 e® 


a” 


4 


+ 


A, k-1 eid => 0. (29) 


It is seen that it is possible to go over from the relations with the upper 
signs to the relations with the lower signs by means of the substitution 


Ark — (— 1) "A_+k 


$ ae es (30) 
Hence, if one has obtained the formula for a A;,, the same transfor- 
mation allows the immediate calculation of A_,+ x. 


From the relations (26)-(29) we have calculated the general expres- . 
sions for the first three coefficients. We refer to the appendix for an example 
of the detailed computation of one of the coefficients. 


Aym-+2,0 = e-t(m+1) 3 (31 a) 
Asmu,0 = e-ims (31 5) 
Asmi2,1 = me et +. 5m + 1) (2m + 3)etim+ys (3) c) 
Qs aa , ip y 
Asmu,1 = — 2(m /‘ 1° tems) 8 4 13 (8m® + 12m + 3) e-ims 
m ey 
sy eas (31 d) 
, 5 : 
Aani,2 = — dei (M41) d __ r Gn ED) ent (m+1)3 
m(m + 1) 33 mas 
+ 96a," e-i (m-1) 


+ 5g, (m + 1) (4m* + 10m + 5) e-ims 


m= - (m+1) (2m+3) (2m+5) (5m-+4) et msya 
(31 e) 
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9 
’ ao” : m(m—1)_; 
A. . Ze ee e-ims ee e-!ms ws . e-tim—2)8 
amt, 2 8 4(m + 1) ‘  480a,? 
__ ipa, 8m? + 12m + 7 ei ms1d 
24 m+ 1 


| m (8m? + 12m + 1) et (m-1) 8 


es . 
1440, 


— 4 (80m! -++- 304m +- 328m? -+- 113m 4-15) e-ims, 


(31 f) 


From those formule we may obtain the expressions for the intensities 
of the spectra, as far as the third term included of the series expansion. 


The final formule are: 


Teemt1 = a {! + q- aa it ba.) cos ¢ 
+ @[-4—F Genet tala 
_ ‘se, x [) 28 25 + a cos 2é 
T ; ra iT 130, ain é 
—- a (256m? + 352m? + 92m + 15) WY ediinciate } P 
(32) 


~~ gers ome m ‘ 1 
Tiamte ee J2m+2 [((m +1)IP {1 <b é —_ cos 8 


aie 


4. ¢2 =e See a, ¥ (m + 1)? 

+] 3 — aoa + Tea 

+s ) cos 25 + 1a," + 1)sin & 

— ffs (m + 1) (2m + 3) (8m + 5)| + esi . (33) 


Those series show an asymmetry in the spectra with regard to the zero 
order. The asymmetry is contained in the second term and is not affected 
by the introduction of terms in p; hence, there is consistency with the results 
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of the preliminary theory, in which for p=0O and even values of n, an 
asymmetric spectrum is obtained. It is also seen that the symmetry is 
re-established when the phase angle 5 = (k + 4)z, confirming the result 
of Murty and the conclusion of Section 4. 

Although the series (32) and (33) are cut off after the third terms, limit- 
ing their practical use to values of é < 1, it is possible to obtain some pro- 
perties from them, which agree qualitatively with the experiment. For 
m= 0 they show that the intensity of the order — 1 is stronger than the 
intensity of the order + 1, and that the intensity of the order + 2 is stronger 
than the intensity of the order — 2, facts that are confirmed qualitatively 
by the experiments of Ramachandra Rao.® 

Taking for instance a, = 1, the same features may be established for 
m=1; for m=2, however, J; = I, (in the considered approximation) 
and from then on all the intensities of positive orders should be stronger 
than the intensities of corresponding negative orders. 


It has already been remarked by Ramachandra Rao that the nature of 
the asymmetry in the case of parallel sound waves for even values of n is 
entirely different from the asymmetry caused by the oblique incidence of 
the light, in which the intensity of any positive order is greater than the 
intensity of the corresponding negative order. 


7. Expiicir CALCULATION OF THE INTENSITIES BY THE SERIES 
METHOD IN THE CASE n = 3 


Taking into account the form of the series (23), we put for n = 3, 


Po — z Agké* 


k=0 


dismti1 = Sea >} Assmti, ké* 
(34) 


EM+2q, m+1 


di smte ~ 2m+2(m + 1)! 1)! > Axsmsa, ke” 








EM+1y M41 ke 
d+sm+s ial 2m+1 (m + 1)! y | Assmt+s, k& 
k=0 J 


The factors preceding the series “have again been chosen as the corres- 
ponding lowest power coefficients of the series expansion of the solution 
(17) in the case p = 0, excepted for the factor of $43m+2, which should read 
m “+1 1 _ eH *) bt EM+2g +1 


Ge 23 (m + 1)’ 
A6 
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but which has been chosen in the above form in order to simplify the rela- 
tions between the series coefficients. Substitution of the series (34) into 
the equation (9) leads to the following relations between the coefficients 
Ark; 





(m +k + 2) Assmii,ku + Atom, ky + Tmr Axsmia,k-1 
F MAsamx2, ket? + Teen Assmsa,k-1et” 
=23m+1)*Assmsyk  (m> 1) (35) 
(k + 2) Assen F ote + Aasea FS Age naet 


2 ; 
+ + Ax, ket? 
=F Asy (36) 
m+ 1 
ag 


(m+ k + 3) Assmse,ku + 





Assmti, ku + Atsm+s, kit 
ag? 

F (m + 1) Agsmzi1€F?% + 4(m + 2) Azsmss, net" 

=F Gm+2)*Arsmiek (m>1) (37) 


1 . 
(kK + 3) Asoku + q, 421, ku Axs, ku + Axgikuer”® 


a3” +18 
+ 7 Ais, k-1€ 


= 2ipA+2,k (38) 


(m+k + 2) Asemas, kt + 4 Assma, kat Aism+a, k-1 


, 2 : 
F (m+ 1) Azam, ky eF + 4 mid) Aismse, n-e*"” 


= 5 (Gm + 3)*Assms,k — (m>0) (39) 
(k + 1) Ao, kH — t A_;, k-1 + 4 A,, k-1 
_ os" 


= 2 P 
4 A_s,K-1e~% + 7. As, k-1 e3 = 0, (40) 





= 5 





ET 











sla- 
nto 
nts 


—1 


35) 


36) 


+1 


£18 


38) 


$18 


39) 
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In order to transform the formule with the upper signs into the formule 


with the lower signs the substitution: 


Ark =“ (— 1)"A+,k \ 
8 > —$§ 


From the relations (35)-(40) one obtains, 





(41) 


must be performed. The same transformation changes the formule for 
A,~ into the expressions for A_,, x. 


Agmis, g = e-t (M+1) 8 (42 a) 
Asm, gq = e-ims (42 b) 
1 . 

Asmi2,0 = (AK ~ ci?) ied (42 c) 
Agmis,1 = ¢ip (m + 1) (2m + 3) ent mins (42 d) 
Asmea,1 = $m + 1) (6m + etm (42 e) 
Asmia,1 = 52 (m + 1) (18m? + 33m + 10) ei 

a3 

— © (6m + Lm + 6) etemenra (42f) 

+t oe 1 ee a,” -i 
Agmis, 2 = Ya, e 4° + (m+1) 4(m +2)° (M+1) 6 
9 
— Teo?” + 1) (2m + 3) (2m + 5) 
x (5m + 4) et m+n) 8 (42 g) 
mM giumays_! 4 
Asm, s = 36. ss e i ( 1) 8 e ims + g ws — (m+ 1)° (M+1)8 
a. ae e-ims aa, 
— 4m FD H° — 480 (540m* +- 1692m$ 
+ 1305m? + 213m + 20) em (42 h) 


pr st, (m + 1) (6m® + 15m + 7) e-ima 


+ terms in ipe? «m+ 6 


(42 i) 
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A _ Ipas 18m? + 21m + 13 i miss ; 
3M+1, 3 ~ <a; 0° : 


i. ie), ek eee 


+ Rae lo,” (6m? + 7m — 1) et m-008 


+ terms in ip e-*™, (42 j) 


In (42) and (42 /) only those terms are written that give a contribution 
to the term in é* in the series for the intensities. 


We have finally calculated the expressions for the intensities as far the 
first term of the series contributing to the asymmetry. We find 


E2M+2, 2m 1 aq” 
l¢smti = 2amF8 m1? {1 + #|- 2 2(m +1) 
lfm a3 
+ 4 (ime, + m1) 8? 


— ia (432m + 390m? + 3m + 5) | 





Sas 
m+ 1 
1 E2M+4, .2M+2 1+ (m+1)? m+1 

t3mt2 = 79m+4) [(m + Dip 4a,” ag 


Z » ests 


ag 


taf + sind + ....} (43) 





cos 5 


F €3 





sind +....} (44) 


EeM+iq, 2M+2 


& 
lismt3 = 22m+2 [m+ 1 ppt va 5[! + a5 


1 
7 cos 8 + fh p? (m + 1) (2m + 3) 





a # 
x (8m + 5)|+€ aaa, (m + 1) sin 8 + oe : 
(45) 
The intensities of those spectra also show an asymmetry with respect | 
to the zero order line. The asymmetry disappears when p is put equal to 


zero, thus confirming the results of the preliminary theory, or when the 
phase angle 6 = km in accordance with the conclusions of Section 4. 


From the formule it is to be expected that at least for values of € < 1, 
the intensities of the orders 3m-+ 1 and 3m-+3 should respectively be 





on 


the 


43) 











q 
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greater than the intensities of the orders — 3m—1 and — 3m — 3, and 
that the intensity of the order — 3m — 2 should be stronger than the intensity 
of the order 3m-+ 2. Here also we have a marked difference with the 
nature of the asymmetry in the case of oblique incidence of the light. 


With regard to the lowest power term in é we may also notice that for 
small values of é the intensities of orders + 3m + 2 should be much weaker 
than the intensities of the orders + 3m+1 and + 3m+3. 


We wish to express our sincere thanks to Prof. Dr. M. Nuyens for his 
continuous interest in our work. 


SUMMARY 


The theory of the diffraction of light by two superposed parallel 
supersonic waves, consisting of the fundamental tone and the n-th harmonic 
is developed, starting from the wave equation for the electric field of the 
light. The Fourier series method, first used in Raman and Nath’s general- 
ized theory, is here employed to derive a system of difference-differential 
equations for the amplitudes of the diffracted light waves. The r-th order 
spectrum makes an angle 0, = — Arc sinr A/A* with the direction of the 
incident light and presents a change of frequency — rv*. In the case that 
the right-hand side of the difference-differential equations may be neglected, 
the exact solution is obtained by means of a complex function method. 
From the structure of the general system of difference-differential equa- 
tions it is shown that the intensities of the orders r and —r are always dif- 
ferent for an even as well as for an odd ratio of the sound frequencies, 
excepted for some special values of the phase angle of the sound waves. 
A solution of the general system for m = 2 is written in the form of a power 
series in ¢, the terms of which are calculated till the third ones included; 
the asymmetry of the intensities of opposite orders is not due to the terms 
in p. In the case n = 3 a series solution also leads to an asymmetric pattern 
with respect to the zero order line; the terms in p are here responsible for 
the asymmetry, so that the symmetry property reappears for p=0, in 
accordance with the simplified theory based on Raman and Nath’s preli- 
minary theory. 
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APPENDIX 


Example of the calculation of a series coefficient from the relations existing 
between them: 


Firstly we must remark that the boundary condition ¢, (0) = 0 always 
requires Aggy =1. The second boundary condition ¢,(0)=O(r40) is 
always satisfied, a power of € greater than zero always preceding the series 
as a factor. 


We shall now give a typical example of the calculation of a coefficient 
in the series ¢?ami2, namely Agmy2,; The principle is quite the same for 
the establishment of the other coefficients, those for n = 3 included. Only 
the length of the computations increases considerably as one tries to obtain 
higher coefficients. 


Putting into the relation (28)k =0, observing that the coefficients, 
the second indices of which are negative, should be identically zero, and 
taking into account the values of Agms;,9 and Agmis, o (31 a, b), previously 
calculated by the same method, one has 


(m + 2) Aemiz,1 — (m + 1) Aam, so” 


= MEF e-ims + 2 (m + 1% ipetomn, 
Qe 


We write this relation for m=0, 1, 2,...... , m—1, m: 
= 1 25, aid 
2Ao — 1Ay,€ aaah. =* + 2.1 Ipe 
2 


3Ag — 2Age® = e+ 2.2% iperiv 
2 


dey — 3A ye? = > e-i® 4 2.32 iperio? 
Qe 


(m + 1) Agm,1 — MAgm,1e-? = ia, m8 + 2.mtipe™ 
(m + 2) Asmse,1 — (m + 1) Aam, so% 


= me etd + 2 (m + 1)tipet mend, 
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From (29) we get for k = 0, 4g, =0. Multiplying both sides of the 
above equations respectively by e-*™, etim-8, e-tim-2)8, eB ang | 
1, and adding them, we obtain 


m+ 1 m+1 


(m + 2) Agm2, 1 = a, > pe-ims a 2ip ; p? e- (M+1) 3 
p=1 


p=1 


Using the formule 





_n(n+ 1) P _n(n+1)Qn+ 1) 


p=1 pel 


there comes finally 


m+1 enims , (m + eet 3; 


Agm2, : => “4 


ipe- (M+1) 3° 
a2 

It is perhaps interesting to remark, that if one wishes to calculate a 
great number of series coefficients one may first calculate Ajo, Aip,. 
Agpi1, 9 (p being a given number), then calculate Ag;, Au,.-.-, Mep,1, further 


compute App, Ay, ...-, Aep-i,e, etc. Hereby the relations (26) and (28) 
are used alternatively. 


Analogue methods may be used for the direct’ calculation of the 
coefficients in the casen = 3. This method is particularly suited, for given 


values of the parameter p and the phase angle 5, when using a computing 
machine. 
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